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I . Introduction 

The goal of the research is to develop a unified constitutive model 
that is applicable for high temperature superalloys used in modern gas 
turbines. The formulation will be considered successful if: (1) the 
resulting formulation is efficient for numerically intensive computation 
such as found in nonlinear finite element models, (2) there is a direct 
correspondence between the material parameters and experimental data, and 
(3) the resulting formulations are reasonably accurate for a variety of 
loading conditions. 

Two unified inelastic state variable constitutive models have been 
evaluated for use with the damage parameter proposed by Kachanov [1]. The 
first is the model of Bodner and Partom [2,3] in which hardening is modeled 
through the use of a single state variable that is similar to drag stress. 
The other constitutive model is an extension of the Bodner-Partom flow 
proposed by Ramaswamy [4] that employs both a drag stress and back stress. 
This extension has been successful for predicting the tensile, creep, 
fatigue, torsional and nonproportional response of Rene' 80 at several 
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temperatures. In both formulations a cumulative damage parameter is 
introduced to model the changes in material properties due to the formation 
of microcracks and microvoids that ultimately produce a macroscopic crack. 
Calculations have shown that the drag stress/damage model is reasonable for 

predicting the tensile and creep responses of IN100 at 1 350°F and Rene’ 95 

at 1200°F, but the model is not entirely satisfactory for predicting the 
cyclic response of these materials. In this study a back stress/drag 
stress/damage model has been evaluated for Rene' 95 at 1200°F and is shown 
to predict the tensile, creep, and cyclic loading responses reasonably well. 
II . Drag Stress/Damage Model 

The initial phase of this research is based on a constitutive model 
using only one hardening variable, Z, to simulate the drag stress of the 
material and a damage variable, u, to model the changes in the material 
properties due the formation of microcracks and microvoids. The inelastic 
flow equation of Bodner was used in the following form 


■ij 


/ n+1 ,Z*(1-oj)^n, 


S. . 
iJ 
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where is the deviatoric stress tensor and Jj - S^j S^j. The value of n 
controls strain rate sensitivity and D# is the limiting strain rate. The 


evolution equations for the state variables Z and uj were developed using 
the Hemholtz free energy as potential function similar to the work in [7]. 
The results show that the flow law and evolution equations are 
thermodynamically associated and are not free to be derived independently. 

The damage parameter, as proposed by Kachanov, was defined to exist on 
the interval [0,1]. A value of zero represents no damage and a value of one 
for 0 ) denotes complete failure. Calculations for Rene’ 95 and IN100 
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indicated that a. value for the damage state variable of much less than one 
corresponds to failure. 

The values of the material constants and initial value of the state 
variables were determined from experimental data that was obtained in a 
previous study. The experimental work was conducted by the Air Force Wright 
Aeronautical Laboratory, WPAFB, Ohio and Mars-Test Inc., Cincinnati, Ohio 
[5,6]. The correlation of the constitutive model to experimental 

observations of monotonic tensile and creep tests on INI 00 at 1350°F and 

Rene’ 95 at 1200°F was very good. The calculated response for INI 00 is 
shown in Figures 1 and 2. 

Extending the drag stress/damage model to predict fatigue loop 
responses was not completely successful. The model appeared to have trouble 
capturing the shape of the fatigue loop. This shortcoming is possibly due 
to how the state variable Z quantifies dislocation movements for a 
particular type of loading. Drag stress is metallurgically associated with 
the retardation of dislocation movement due to the interaction with 
precipitates which results in dislocation Jogs or looping for example. Back 
stress is associated with the pile-up of dislocations at a barrier, such as 
a grain boundary, which cannot be overcome. During cyclic loading the back 
stress oscillates with the applied load while the drag stress either 
increases monotonically or remains essentially constant. Since the 
dislocation movement was only in one direction for the tensile and creep 
teats, the contributions of the drag stress and back stress to the total 
resistance of dislocation motion could be simulated by one state variable 
monotonically increasing to a steady state value. Although this lumping of 
resistances to dislocation motion worked well for monotonic loadings, it 
appears inadequate for cyclic loadings. 
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III. Drag Stress/Back Stress/Damage Model 


Based on the metallurgical considerations given above and the success 

of the drag stress/back stress model proposed by Ramaswamy for predicting 

the fatigue response of Rene* 80, a drag stress/back stress/damage model is 

proposed. The drag stress scalar Z simulates long term cyclic hardening 

while the back stress tensor Q. . models the short term strain or work 

^ j 

hardening. The governing equations are: 

Inelastic Flow Equation, 


^ij » Do exp 




Back Stress Evolution Equation, 
I 




ij 




Drag Stress Evolution Equation, 
Z - m(Zi-Z)W^ 

Damage Evolution Equation, 


( 2 ) 


(3) 
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and 

Back Stress Relaxation Equation 


“3 ■ <V“33t> 




Where K 
f, are material constants. 


( 6 ) 

). In the back stress evolution equation fj and 
In the drag stress evolution equation the value 


of m controls the hardening rate, Zj is the saturated value of Z, and is 
the inelastic work rate. The parameter g* in the damage evolution equation 
is a material function and damage growth is assumed to be proportional to 
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the maximum principal tensile inelastic strain rate e . Initially damage 

0)€tX 

is assumed to be absent in the material before loading. In the back stress 
relaxation equation is the steady state value of the back stress, is 

the maximum saturated value of the back stress observed in uniaxial tensile 
tests, the constant Oo is introduced to nondimensionalize stress, and B and 
r are material constants to control the time dependence. The proposed model 

reduces to that of Ramaswamy for u » u » 0. 

In uniaxial loading Rene* 95 response displays tension-compression 
asymmetry as shown in Figure 5. At constant strain rate the compressive 
tests saturate at a higher magnitude of stress than do the tensile tests. 
The experimental fatigue response as given in Figure 9 for example shows 
that for the same level of stress the tensile strain is less than the 
compressive strain. It is also shown in [5] that the minimum creep rates 
observed in compression are much less than those in tension at corresponding 
values of stress. These pieces of information are used to substantiate that 
a higher initial value of hardness for the material in compression than in 
tension. For the uniaxial exercise this asymmetry was Included in the model 

•f - 

by having two drag stress state variables, Z and Z similar to kinematic 
hardening components. The corresponding evolution equations are: 

Tensile Drag Stress Evolution Equation, 

Z = m (Z,-Z )W^ (7) 

and 

Compressive Drag Stress Evolution Equation, 

Z~ = m"(z7-Z~)W^. (8) 
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When the stress is greater than zero Z*Z^ and when the stress is less than 
— • + 

zero Z=Z . The Z evolution equation was only active for stresses greater 

• — 

than zero and the Z equation was only active for stresses less than zero. 
IV. Evaluation of Material Parameters 

The values of Do. n, D , Z* and Z, were found using a nonlinear 

nia X 

regression analysis of the inelastic strain rates corresponding to saturated 
stress levels observed in the tensile tests. The values of the material 

constants fi, fj, m^and m were determined using an iterative computer 
program with the tensile, creep, and cyclic test data as inputs. The 
program iterated until it converged on stable values for these material 

constants. For Rene' 95 at 1200®F the values of Z* and Zi are equal. 

The damage function, gj, was determined as the final step after all 
other constants had been found. The damage growth was evaluated from the 
tertiary creep response for small values of accumulated inelastic strain and 
the long term fatigue response corresponding to large values of accumulated 
inelastic strain. Using both the creep and cyclic tests the values of gj, 
relating damage to the accumulated inelastic strain, were found as shown in 
Figure 3* This damage curve demonstrates that the greatest damage growth 
occurs early in the loading of the specimen and the rate of damage 
accumulation approaches a steady-state value rather rapidly. 

The maximum value of the back stress is observed in the short time 
tensile tests when the stress saturates. During long time loading such as 
creep, the maximum value of the back stress has a lower value than that in 
the short time tests. This lower value was seen by Ramaswamy [4] for Rene' 
80 at 1^400“F and 1600*F. This reduced back stress value is also present in 
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Rene' 95 at 1200°F. Figure shows how the maximum back stress decreases 

from the maximum value as a function of the applied stress. Without 

test data to for substantiation a lower bound has been artificially imposed. 
The value of Oo was arbitrarily chosen as 200 Ksi to nondimensional ize 
stress. Since the time required for a given creep test to reach a stable 
back stress can be observed, the values of B and r were determined by a 
linear regression of the back stress relaxation equation and the teat data. 
V. Comparisons of Experimental and Calculated Results 

The strain rate control test data for Specimens 1-1 and 6-1 in Figure 6 
show that the saturation stress for Rene' 95 is essentially strain rate 
independent. The model predictions show a small amount of strain rate 
dependence; however the predictions are reasonable for both the stress rate 
controlled and strain rate controlled tests. 

The comparison of the observed and predicted creep response is shown in 
Figure 7. The creep stresses range from 1H6 Ksi to 175 Ksi. The 
experimental data is ordered except for the 1*16 Ksi test. The model 
predicts reasonably well the primary, secondary, and tertiary creep 
responses for all of the tests except the 1 M6 Ksi test. The slight dip in 
the predicted creep response connecting the secondary and tertiary creep 
regions is a consequence of matching the constants in the back stress 
relaxation equation with the damage growth equation. However, once the back 
stress is stabilized the model predictions are parallel to the experimental 
creep responses. 

A typical stress control cyclic test is shown in Figures 8 and 9. The 
predicted and experimental peak strains as a function of time are shown in 
Figure 9. The experimental tensile peak strain is observed to remain 
relatively constant while the compressive peak strain decreases with each 


193 




additional loop until a steady state response is reached. The observed 
error in the model prediction is due to a combination of factors. However 
since test 5-3» Ref. C5]. is for a stress range of ± 168 Ksl the deformation 
is well into the plastic range. It can be seen from the tensile tests in 
Figure 6 that a small deviation in the peak stress produces a large 
variation in the predicted tensile peak strain. Since the trends are 
clearly correct the predicted tensile strain is considered acceptable. 
Figure 8 shows how the shape of the fatigue loop changes with increasing 
cycles. A displacement limit of 0.02 In/in was defined as failure since the 
model does not include a criteria for the transition from damage 
accumulation to crack growth. The model has been used to predict the 
fatigue life up to cracking for two stress control tests 5-3 and 3-6 of Ref. 
[5] at ±168 Ksi and ±158 Ksi, respectively. In these calculations the 
damage growth appears to be reasonably accurate. Additional fatigue 
calculations will be made in the future to further evaluate the model. 
Finally, future work to use the model to predict the material deformation to 
combined creep and cyclic loadings is planned. 
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Figure 1. Predicted and Experimental Tensile Response of IN100 at 1350®F 
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Figure 2. Predicted and Experimental Creep Response of IN100 at 1350®F 
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Figure 3 
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Damage as a Function of Accumulated Inelastic Tensile Strain for 
Rene* 95 at 1200°F 



Variation of Saturated Back Stress with Applied Stress for 
Rene' 95 at 1200®F 
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Figure 5 


Figure 6 
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. Predicted and Experimental Tensile Response of Rene* 95 at 1200°F 
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Figure 7 


Figure 8 
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Predicted and Experimental Creep Response of Rene* 95 at 1200®F 
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Variation of Fatigue Loop for Specimen No. 5~3 Loaded at ±168 Ksi 


at 10 Cycles per Second for Rene* 95 at 1200°F 
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Figure 9 



Predicted and Experimental Strains for Specimen No. 5-3 Loaded at 
±168 Ksi at 10 Cycles per Second for Rene* 95 at 1200®F 
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